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Abstract 

We proposed an optimal asymmetric resource allocation 
(ARA) scheme for the decode-and-forward (DF) dual-hop 
multi-relay systems in the downlink, with Orthogonal 
Frequency-Division Multiple Access (OFDMA) for multi- 
user transmission. Our work is different in that the time slots 
for the two hops via each of the relays are designed to be 
asymmetric, i.e., with K relays in the cell, a total of 2 K time 
slots may be of different durations, which enhances the 
degree of freedom over the previous work. Also, a 
destination may be served by multiple relays at the same 
time to enhance the transmission diversity. Moreover, 
closed-form results for optimal resource allocation are 
derived, which require only limited amount of feedback 
information. Simulation results show that, thanks to the 
multi-time and multi-relay diversities, the proposed ARA 
scheme can provide a much better performance than the 
scheme with symmetric time allocation, as well as the 
scheme with asymmetric time allocation for a cell composed 
of independent single-relay sub-systems, especially when 
the relays are relatively close to the source. 
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Introduction 

Recently, relay technologies, including amplify-and- 
forward (AF) [Yang08] and decode-and-forward (DF) 
[Kim08] inodes, have been considered as a promising 
and achievable solution for the fourth-generation (4G) 
mobile communication systems such as Long Term 
Evolution Advanced (LTE-Advanced) systems 
[Yang09] [IwamuralO]. 

Adaptive resource allocation (RA) [Shen05] plays an 
essential part in relay systems. In [Li08], subcarrier 


pairing and power allocation was conducted 
successively in a single-relay system for both AF and 
DF modes to maximise the system capacity. In 
[Wang08], the optimal joint subcarrier matching and 
power allocation was proposed for DF relay systems 
with single destination. However, it requires the bits 
to be grouped in relaying, moreover, the complexity of 
subcarrier pairing is too high to guarantee that the 
channel state remains the same during this process. In 
[Nam07] an optimal RA algorithm was proposed to 
maximise the bandwidth efficiency of the DF relay 
system, however, it assumes fixed geometric locations 
of the relays, which is merely practical in the real 
urban environment. While in [SalemlO], the authors 
considered a multi-relay system with random relay 
locations, and also proposed an optimal RA algorithm. 
However, in most previous work, the time slot 
durations for the two hops were designed to be 
symmetric. In [Agustin09], a two-way relay channel 
model was investigated, and an asymmetric time 
allocation scheme was proposed, where the two 
consecutive time slot durations were designed to be 
asymmetric. However, it is only assumed that a 
uniform power allocation is conducted across all 
subcarriers, which was not efficient. In [Zhoull], an 
optimal ARA scheme was proposed for an OFDM 
based DF relay system, however, it can only be 
applied to a scenario where a cell was divided into 
multiple single-relay sub-systems with RA conducted 
in each of individual sub-systems independently. 
Hence, the ARA scheme is not readily extendable to a 
multi-relay multi-destination system. 

In this paper, we investigate an ARA scheme for DF 
dual-hop multi-relay OFDMA systems. Our work is 
different in the following aspects. Firstly, to the best of 
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our knowledge, this is the first work to apply 
asymmetric time allocation to a general multi-relay 
multi-destination system, where the time slots for the 
two hops via each of the relays are designed to be 
asymmetric, i.e., with K relays in the cell, a total of 2 K 
time slots may be of different durations. As a result, it 
enhances the degree of freedom for transmission over 
the previous symmetric RA (SRA) scheme. Secondly, 
an optimal algorithm is proposed to perform joint time, 
power and subcarrier allocation to obtain the global 
optimal results, with only limited amount of feedback 
information from relays and destinations. This is not 
an easy extension of [Zhoull] where each destination 
is served by only one relay in an independent sub- 
system. While the proposed work allows multiple 
relays to serve a single destination, which enhances 
the degrees of freedom. Simulation results show that, 
thanks to the multi-time and multi-relay diversities, 
the proposed ARA scheme outperforms the SRA 
algorithm in [SalemlO], as well as the ARA algorithm 
in [Zhoull], with higher achievable system 
throughput, especially when the relays are relatively 
close to the source. Moreover, we demonstrate the 
impact of the relays' locations on the results of 
asymmetric time allocation. 

Section II presents the system model and the problem 
formulation. In Section III, the optimal ARA algorithm 
is proposed. Simulation results are shown in Section 
IV, and the conclusion is drawn in Section V. 

System Model and Problem Formulation 

System Model 

We consider a DF OFDMA cellular system in the 
downlink, with two hops and K relays, as illustrated in 
FIG. 1, where S denotes the source node (BS), 
D ( (/ = 1,..., Lj denotes the different destination nodes 

(mobile users), and R k (k = denotes the 

multiple relay nodes. The K relays are uniformly 
distributed within the range of width w in FIG. 1. We 
define d SR as the distance from the source to the mid 

circle of relays' band. Without the loss of generality, it 
is assumed that all the destinations have the 
approximately same distance from the source, which 
is denoted by d SD , and they are uniformly distributed 

along a circle. The hop from source to relays is defined 
as the first hop, and the hop from the relays to 
destinations is the second hop. We define { I , . . TV } as 
the orthogonal subcarriers set, where N is the total 


number of subcarriers. 
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FIG. 1 DUAL-HOP MULTI-RELAY CELLULAR SYSTEM ( S - 
SOURCE; R - THE k TH RELAY ( k = 1, ..., K ); D - THE l TH 
DESTINATION ( l = 1, ..., L )) 

The total transmits time duration is denoted as T. The 
time slot durations via the relay k for the first and the 
second hops are T k , and T, k , respectively. Assuming 
perfect channel estimation at both relay nodes and 
destinations, we define h!f k as the normalised channel 
impulse response on subcarrier n between the sources 
and relay k in the first hop, and h^l , as the normalised 
channel impulse response between relay k and 
destination l in the second hop. The path loss (PL) 
attenuation coefficients in the first and the second 
hops are denoted as Aj and A, , respectively. We 

define d k as the distance between the sources and 
relay k in the first hop, and d kl as the distance 

between the relay k to the destination / in the second 
hop. The PL exponents in the first and second hops are 
denoted as r 1 and r 2 , respectively. We define pf ] k as 
the transmit power for the source to relay k on 
subcarrier n in the first hop, and p^l , as the transmit 
power between relay k and destination / in the second 
hop. e {0,l} denotes subcarrier allocation 

indicator for subcarrier n in the first hop. If subcarrier 
n is allocated to relay k, | 1 < " ) = 1 , otherwise, = 0 . 

Similarly, , is the subcarrier allocation indicator in 
the second hop. 

The channel-to-noise ratio (CNR) is 

Yi'k = ^|^i"(r ) | /(AW 1 ) i n the first hop, where N 0 

is the single-sided power spectral density of the 
additive white Gaussian noise (AWGN). The 
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bandwidth for each subcarrier is B/N . In the second 


as the 


hop, we define = n|^ ( |Y(A 2 N 0 B^) 

CNR. The maximum achievable throughput (in bits 
per second) for relay k in the first hop can be written as: 


Q,*=^Z£>g 2 ( 1 + pM) 


( 1 ) 


and the maximum achievable throughput of the 
second hop can be written as: 


■n L N 

^=^ZZ4V°g 2 (l + P^’,) 

IV 1=1 o= 1 


( 2 ) 


Hence, the overall maximum achievable throughput 
during time slot T on the route via relay k can be 
written as: 

TC k = min {T u C u ,T 2#i C 2 J (3) 


£c =i ' zze=^ m («) 

*:=i *:=i (=i 

Constraint (5) implies that the time slot durations for 
both hops via each relay may be asymmetric, while in 
[SalemlO], they are assumed to be equal, i.e., 
Tik = T 2k = T/ 2V(n) . Constraint (8) indicates the 
exclusivity of subcarriers allocated to nodes. 

Optimal Asymmetric Resource Allocation 

In this section, we propose the optimal solution to the 
objective function (4). It can be derived that (4) is 
optimised when T lk C lk =T 2k C 2k (yk) [Zhoull]. Also 
using (3), (4) can be rewritten as: 


k 1 

max V— (r. .C. 

v'$A n L.\k^ t fa 2 v ' ' 


k + T 2 ,k 


C 2,k) 


(9) 


Problem Formulation 

In this subsection, the RA problem of the investigated 
system is formulated. We define P T as the total 

transmit power of the system in the downlink, P s as 
the transmit power of the source node, and P R as the 


It can be shown in the Appendix that (9) is a convex 
function. Thus (9) can be solved using the Karush- 
Kuhn-Tucker (KKT) conditions [Boyd04]. Using (9), 
and (5)-(8), the Lagrangian function can be built as 
equation: 


total transmit power of all the relay nodes. We also 
assume that each subcarrier is occupied by only one 
node in a time slot. Unlike SRA in [SalemlO], where 
the time slot durations for the two hops are identical, 
the proposed ARA scheme allows all the 2 K time slots 
to be of different durations. Also, it allows a certain 
subcarrier to be allocated to any of the K relays in the 
first hop, while in [Boyd04], N subcarriers are divided 
into K groups and each group serves a single-relay 
sub-system. In the second hop, our system allows 
multiple relays to serve one destination, which is 
different from the multiple single-relay sub-systems, 
where each destination can only be served by the 
single relay. We maximise the achievable system 
throughput, which is formulated by: 

K 


K 1 

k = 1 z 


T u EC lo g 2 ( 1 + PuV$) 

n=l 

^ZZC log^/^) 


-Z4 


N 

^uZC lo g 2 ( 1 + PuVu) 


L IV 

- T 2,*£X4w l0 g2( 1 + PMu) 

1 = 1 n = 1 

( K N \ f K L N 

X2>u- p s -aJ 


k = 1 1=1 n = 1 

f K 


-Z a k ( T u +T 2,k-T)-H s i n) Z C - 1 

k = 1 

-z^fzz^-i 


n= 1 V k=l 


(10) 


max Yj C k 

(4) 

where constant B/N is omitted in order to simplify 

„(«) „(») T T ““ 

P\,k ’P2,k,l' 1 \,k' 1 2,k k= 1 


subject to: 


the derivation. We define (fc = 1,...,K) , /./, , //, , 

T u+T 2 ,=T, (Vfc) 

K N K L N 

(5) 

o’,, (^ = 1,...,X), <J 1 ( " ) , and J 2 " ) (n = l,...,Al) as Lagrange 
multipliers. 

ZZpu^ p s' ZZZpS,^^ 

/c=l n=l A:=l Z=1 n=l 

(6) 

Each subcarrier is allocated to one node exclusively. In 
the first hop, the subcarrier h is allocated to relay node 

where p\ n) k > 0, p^h- '0, (V M,n ) 

(7) 

k' when k* satisfies the following rule, which can be 

P +P =P 

S R T 

derived from the KKT condition [Nam07], as: 
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k = arg max 

k 

M 


— /I, 


i°g 2 (i +P%r%)-W% 


(ii) 


In the second hop, the subcarrier n is allocated to the 
relay-destination pair |fc , 1 j which satisfies the same 
rule, as: 

k , l | = arg max 

M 

- Wzli]'{v k ' 1 ) 


-+A 


log 2 (l + M) 


( 12 ) 


Note that k and k may denote different relays. 
Assuming A , , A . , //, , and //, are the optimal 

Lagrange multipliers, the resulting optimal power and 
time slot durations allocated to relay k and k can be 
derived as: 


p (n \ 
r 1 ,k 


(iiVl 

/j 1 In 2 


y<"> 

'l,Jt 


P (H) . = 
r 2 ,fc ; 


(tHf i 

/t In 2 y (n) . 

2 ' 2 ,k ,1 


T , = 

1 ,k 


T . = 

2 ,k 


— K- 

2 k 


- + A, 
2 k 


(13) 

(14) 

(15) 

(16) 


where [*T = max|*,0} . Note that when symmetric 

time allocation [SalemlO] is applied, i.e., 
T . = T . = T/2 , the optimal power allocation results 

are the same as (13)-(16), except that ^l/2-A . j in (13) 

and (l/2 + A. j in (14) are replaced with (l/4 - /L j 

and (l/4 + /L j , respectively. Hence, the proposed 

ARA algorithm has a similar complexity to the SRA 
algorithm. In addition, the value of these multipliers is 
the limited information which needs to be fed back 
from the relay nodes to BS. The amount of feedback 
information increases linearly with the number of 
relays in the cell. 

The following (X + 2) equations are to be solved to 
obtain the optimal Lagrangian multipliers: 


K N 

'i-ZM 

k = 1 n=l 




k = 1 1=1 n = 1 


:(») 


1 ,k 

//j In 2 

s(») 

1 

+ 

>2 ,k,l 

ju 2 In 2 


^15 


= 0 


y(") 

/ 2-Jfc-Z 


= 0 


(17) 


(18) 


( 1 + H)I ZC 


lo g 2 


(i + A.) 2 


h 2 


In 2 


y<"> 

/ 2.Jk,Z 




io g 2 


//j In 2 


y(") 
I l,k 


(19) 


= 0,(Vfc) 


Various search algorithms such as the sub-gradient 
algorithm [Nam07], can be applied to solve (17)-(19). 


Simulation Results 

In our simulations, we use N= 256 subcarriers, K=A 
relays, and L=8 destinations (except for FIGs. 3 and 4). 
All relays are uniformly distributed within a band 
region of width w = 20 m, as shown in FIG. 1. The 
channel is modelled as six independent Rayleigh 
fading paths with the root-mean square (RMS) delay 
spread of 0.5 // s. We apply the Type-D (Roof-to-Roof) 
PL model [Hart06] in the first hop with PL attenuation 
coefficient A 1 = 2.05_f 2 ' 6 x 1 0 26 , where f c = 5 GHz is 
the central carrier frequency. The PL exponent for the 
first hop is r, = 4.5 . Applying the Type-E (Roof-to- 
Ground) PL model [Hart06] in the second hop, we 
have the PL attenuation coefficient as A 2 = 38.4 dB, 

and PL exponent as r 2 = 3.5 , respectively. The total 

bandwidth is B=50 MHz for the downlink [Seidel08], 
and the single-sided power spectral density of AWGN 
is N 0 = -173.8 dBm/Hz. The total transmit power is 

P T = 40 dBm, and the total powers of P s and P R are 38 

dBm. The total transmit time is T= 5 ms. In all Figures 
except FIG. 5, the distance between source and 
destinations is fixed at d SD = 2000 m. The distance 

between the source and the mid circle of the relays d SR 

varies between 600 m and 1800 m. The ARA approach 
in [Zhoull] is implemented for comparison, where the 
whole cellular system is divided into K= 4 single-relay 
sub-systems. 

FIG. 2 illustrates the impact of the locations of relays 
on the maximum achievable system throughput. The 
proposed ARA scheme outperforms the single-relay 
based ARA scheme in [Zhoull] and the multi-relay 
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based SRA scheme [SalemlO], especially when the 
relays are relatively close to the source. For instance, 
when d SR = 600 m, the performance of the proposed 
ARA is around 46.9% and 70.5% higher than that of 
single-relay ARA and SRA, respectively. It can also be 
observed that the performances of proposed ARA and 
single-relay ARA have the same tendency of 
increasing when d SR is in the range of 600 m to 1800 m. 



FIG. 2 IMPACT OF DISTANCE BETWEEN THE SOURCE AND 
THE MID CIRCLE OF RELAYS ON MAXIMUM ACHIEVABLE 
THROUGHPUT, WITH K = 4 RELAYS AND L = 8 DESTINATIONS 

When the d SR further increases, the performance of 
SRA deteriorates rapidly, as the symmetric time 
allocation cannot compensate for the performance 
imbalance caused by different levels of PL attenuation 
in the two hops. 



FIG. 3 IMPACT OF THE NUMBER OF RELAYS ON MAXIMUM 
ACHIEVABLE THROUGHPUT, WITH L = 16 DESTINATIONS 

In FIG. 3, the impact of the number of relays on the 
maximum achievable throughput is illustrated, where 
d SR = 1200 m. The number of destinations is fixed as 
L=16. It is shown that when the number of relays 
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increases from 1 to 16, the system throughput of 
proposed ARA is enlarged by 65.6% due to multi-relay 
diversity, while the throughput of the single-relay 
based ARA scheme [Zhoull] decreases by 9.6% due to 
lack of multi-relay diversity. When the number of 
relays is 1, the proposed ARA scheme reduces to the 
method in [Zhoull]. 



FIG. 4 IMPACT OF THE NUMBER OF DESTINATIONS ON 
MAXIMUM ACHIEVABLE THROUGHPUT, WITH K = 4 RELAYS 

With the same configuration as FIG. 3, FIG. 4 
illustrates the impact of the number of destinations on 
the maximum achievable system throughput, with K = 4 
relays. It can be observed that when the number of 
destinations increases from 8 to 24, the system 
throughput has little increase, which demonstrates 
that the multi-user diversity has less impact on 
performance compared with multi-time diversity and 
multi-relay diversity, as shown in FIGs. 2 and 3. 



FIG. 4 IMPACT OF NORMALISED DISTANCE BETWEEN 
SOURCE AND DESTINATIONS ON THE AVERAGE TIME 
DIFFERENCE BETWEEN TWO SLOTS, WITH K = 4 RELAYS AND 
L = 8 DESTINATIONS 

Defining d SR /d SD as the ratio of the distance between 
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source and the mid circle of the relays band to the total 
distance between source and destinations, we 
demonstrate the impact of d SR /d SD on A T m in FIG. 5, 

where AT ave = XL £ {( T u -T 2k )/T)/K denotes the 

average normalised time difference between two hops. 
It can be observed that when relays are closer to 
source than destinations, the optimal time duration in 
the first hop is smaller than that in the second hop. For 
instance, when d SR /d SD = 0.3 , AT m =-0.78 for 

d gD = 2000 m, which indicates the time duration 
allocated in the second hop can combat the higher 
level of PL attenuation caused by longer distance to 
achieve the same throughputs with the first hop. It is 
also observed that the longer the distances between 
source and destinations are, the steeper the curve is. 
Therefore, when destinations are further away from 
source, asymmetric time allocation plays a more 
important role in system performance. It can be 
deduced that symmetric time allocation may be a 
preferable alternative at d SR /d SD = 0.65 ~ 0.7 

(irrespective of the value of d SD ), due to easier 
synchronisation. 

Conclusions 

In this paper, we have investigated ARA for DF dual- 
hop multi-relay OFDMA cellular systems in the 
downlink. Our proposed ARA scheme performs joint 
time, power and subcarrier allocation to achieve an 
enhanced degree of freedom from the multi-time, 
multi-relay, and multi-user diversities. Simulation 
results show that asymmetric time allocation plays an 
important role in system performance. As a result, the 
proposed ARA scheme can outperform the SRA 
scheme [SalemlO] significantly due to multi-time 
diversity, at a similar complexity, especially when 
relays are relatively close to source. Also, when the 
distances between source and destinations are 
relatively large, the impact of asymmetric allocation is 
more significant. Furthermore, the proposed optimal 
multi-relay ARA scheme outperforms the single-relay 
based ARA scheme in [Zhoull], thanks to the multi- 
relay diversity. Simulation results also show that the 
multi-user diversity plays a less significant role than 
the multi-relay and multi-time diversities in the 
system performance. 

Appendix - Proof of Convexity of the 
Objective Function (9) 

It can be proved that the objective function (9) is a 
convex optimisation problem, by proving the 


convexity of each TC k in (9). For the purpose of 
simplicity in the proof, we define j? = p\ n) k j ^" k and 

^w=P 2 xi/^ 2 xi ' where ^ and ^ are relaxed to 
continuous variables which are in the interval of [0,1]. 
Using (1), (2), (5) and T lk C lk = T 2k C 2k , we can derive 


T u C u XmXL&M 1 + Q&) 

T V C u XXC lo M 1 + €^) 

And using (1), (2), (20) and TC k = (T lk C lk + T 2k C 2k )/2 , 


we have 


C„=- 


1 c. 


2 { C l,k +C 2 ,k 

1 1 

c c 

^2,k 


c +c 

^l,k ^2 ,k 


Hence, the maximum achievable system throughput 
via relay k can be expressed as 




N [lZAPo S2 (i + Cr^) 


The concavity of lo §2 ( 1 + 9uVu) and 

X w X 1, lo §2 0 + ) are proved first. Note 

that XXC log^l + ^y^) can be decomposed into 
N components, and the nth component is 
l°g 2 (l + ) ■ Defining a function 

/(C'PuO = C lo S 2 ( 1 + PM/C) ' the Hessian 

matrix of f^[ n k ,p[ n) k J can be written as 


H (/) = 


d 2 f d 2 f 

d 2 f d 2 f 

¥^Kk d(p^ k ) 2 

(ri:. 1 )’ Pt! 

l (c+pS)' 

l2 M M 

ri,k ,k 

i ^+ p ' Cl ) 2 (€ ) +* > ) 2 
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Defining r = [r 1 ,r, ] as a non-zero real valued number 
vector, we have 


rH(/)r r 


1 

Pu 

1*2 (C+Pu') 2 

[wr 




J 


(24) 


It is obvious that when 0 < < 1 and p[ n) k > 0 , 

rH(/)r T < 0 always holds, which means that H( f ) is a 
negative definite matrix. Thus, it can be concluded 
that is concave. Accordingly, 

X^^log 2 (l + PM) is concave. Similarly, 

ZwZli^"w lo g 2 ( 1 + Pu/w) is concave. Therefore, 
(22) is convex, and (9) is convex objective function, 
which means that the globe optimal solution to (9) can 
be obtained, as given by (13)-(16). 
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